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Self-induced flavor conversions of supernova (SN) neutrinos can strongly modify the flavor depen- 
dent fluxes. We perform a linearized flavor stability analysis with accretion-phase matter profiles of 
a 15 Mq spherically symmetric model and corresponding neutrino fluxes. We use reaUstic energy 
and angle distributions, the latter deviating strongly from quasi-isotropic emission, thus accounting 
for both multi-angle and multi-energy effects. For our matter and neutrino density proflle we always 
find stable conditions: flavor conversions are limited to the usual Mikheyev-Smirnov-Wolfenstein ef- 
fect. In this case one may distinguish the neutrino mass hierarchy in a SN neutrino signal if the 
mixing angle 6is is as large as suggested by recent experiments. 



PACS numbers: 97.60.Bw, 14.60.Pq 

Introduction. — The huge neutrino fluxes emitted by 
core-collapse supernovae (SNe) are key to the explosion 
dynamics and nucleosynthesis [T] and detecting a high- 
statistics "neutrino light curve" from the next nearby SN 
is a major goal for neutrino astronomy [2]. Besides prob- 
ing the core-collapse phenomenon in unprecedented de- 
tail, one may detect signatures of flavor oscillations and 
extract information on neutrino mixing parameters OS]. 

The refractive effect caused by matter [5] suppresses 
flavor oscillations until neutrinos pass through the 
Mikheyev-Smirnov-Wolfenstein (MSW) region in the col- 
lapsing star's envelope [HI [7]. However, neutrino-neutrino 
interactions, through a flavor off-diagonal refractive in- 
dex [HI H] , can trigger self-induced flavor conversions [TOj- 
[T^ . This collective effect usually occurs between the neu- 
trino sphere and the MSW region and can strongly mod- 
ify neutrino spectra |13H15j . although this would never 
seem to help explode the star [T5]. Actually, in low- 
mass SNe (not studied here) the density falls off so fast 
that MSW can occur first, leading to novel effects on the 
prompt burst |17j . 

Collective oscillations at first seemed unaffected by 
matter because its influence does not depend on neu- 
trino energies |13j . However, depending on emission 
angle, neutrinos accrue different matter-induced flavor- 
dependent phases until they reach a given radius. This 
"multi-angle matter effect" can suppress self-induced fla- 
vor conversion [TH] . Based on schematic flux spectra, this 
was numerically confirmed for accretion-phase SN models 
where the density near the core is large [19] . This epoch, 
before the delayed explosion finally takes off, is when the 
neutrino luminosity and the difference between the P,, 
and Dfj^^T fluxes are largest. If self-induced flavor conver- 
sion did not occur and the mixing angle ^13 was not very 
small [20], the accretion phase would provide a plausible 
opportunity to determine the mass hierarchy [31 119j . 

Numerical multi-angle simulations of collective oscilla- 
tions are very demanding [21|. but it is much easier to 



study if such oscillations are suppressed for given density 
profile and neutrino distributions. Self-induced conver- 
sion requires that part of the spectrum is prepared in one 
fiavor, the rest in another. The collective mode consists of 
pendulum-like flavor exchange between these parts with- 
out changing the overall flavor content [TUl [22] • The in- 
evitable starting point is a flavor instability of the neu- 
trino distribution caused by neutrino-neutrino refraction. 
An exponentially growing mode can be detected with a 
linearized analysis of the evolution equations [121 [23]. We 
here apply this method to a numerical accretion-phase 
SN model, for the first time using both realistic neutrino 
energy spectra and angular distributions. 

Numerical SN model. — Our spherically symmetric sim- 
ulation of an accretion-phase SN model was performed 
with the Prometheus- Vertex code as in Ref. 24J, now 
with a 15 Mq progenitor J25j. The transport module 
computes the energy and angle distributions of ly and 
D for all fiavors by a tangent-ray discretization of the 
Boltzmann transport equation |26j . We used 21 nearly 
geometrically spaced energy bins up to 380 MeV and 672 
tangent rays. We do not artificially trigger an explosion, 
but otherwise our model is comparable to Ref. [19J. We 
use several snapshots and illustrate our findings with one 
at 280 ms post bounce. The fiux spectra (Fig. [T]) show a 

excess from deleptonization and a 9^ flux almost twice 
that of i>x (representing any of z^^.r or i^^,t) and average 
i^e, Vf. and Vx energies of 15.3, 18.1, and 16.9 MeV. 

We study neutrino propagation in the free-streaming 
limit, so we can describe the angular distribution by the 
angle da relative to the radial direction at a chosen inner- 
boundary radius R. Actually it is more convenient to use 
u = siv?'dR = (1 — cos^ -dr) / E? ^ which is uniformly 
distributed on < u < 1 if emission is isotropic at a 
"neutrino sphere" with radius R [13 US]. We choose 
R — 44.7 km and show the corresponding u distribution 
in Fig. [Tj Isotropic emission from a neutrino sphere is 
not a good description because neutrinos emerge from a 
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FIG. 1: Flux spectra for our 280 ms SN model. The angle 
variable < m < 1 is based on R = 44.7 km. 

thick layer. The and Vx intensities are similar in the 
radial direction: the excess De flux largely arises from 
its broader angular distribution (larger emission region) . 
Flavor osciUations depend on the difference of the e and 
X distributions, which is small in the radial direction 
(Fig. [T]) . The angular distributions do not cross, although 
in principle there could have been a forward excess. 

In the context of neutrino oscillations, uj = Arn^ /2E 
is a preferred energy variable, where Am^ = (50 meV)^ 
is the "atmospheric" neutrino mass-squared difference 
relevant for 1-3 oscillations studied here. Moreover, 
treating anti-neutrinos formally as negative-energy neu- 
trinos with negative occupation numbers vastly simpli- 
fies the formalism. Flavor oscillations can exchange 
with i/j., leaving the overall neutrino flux unchanged, so 
only F^^—F^^ matters. Our sign convention means that 
for anti-neutrinos we then use Fp^ — Fp^ , corresponding 
to the flavor isospin convention [13j . The neutrino flux 
difference distribution g{oj,u) thus defined is shown in 
Fig. [2] It is negative for anti-neutrinos (w < 0) because 
F[,^ > Fjy^. For u! ~ 0.2 km~^ there is a spectral crossing 
as a function of u, i.e. for large E the i^x flux does exceed 
the Ve flux in the forward direction. 

Self-induced oscillations exchange the positive and neg- 
ative parts of g{u!, u), leaving fixed the overall flavor con- 
tent e = (Fi,^-i^i.J/(i^i>^-FpJ-l = f du}dug{u},u). Our 
g{i^, u) is mostly negative for P and mostly positive for z/, 
so collective oscillations largely correspond to pair con- 
versions Ve^e ^ Vxi^x- Accretion-phase distributions are 
"single crossed" in this sense, i.e. g{LiJ, u) changes sign 




FIG. 2: Distribution g{uj, u) describing the neutrino fluxes. 

essentially only on the line cj = 0, because of the large 
excess of the and i>e fluxes. Significant multiple cross- 
ings are typical for the cooling phase 

Equations of motion (EoM). — We describe three-flavor 
neutrino propagation by energy- and angle-dependent 
3x3 matrices ^E,u{f')- Boldface characters denote ma- 
trices in flavor space. The diagonal ^e,u elements are 
the ordinary number fluxes F^ ^ (flavor a) integrated 
over a sphere of radius r, with negative E and nega- 
tive number fluxes for anti-neutrinos. The off-diagonal 
elements, which are initially zero, represent phase in- 
formation caused by flavor oscillations. The flavor evo- 
lution is then provided by the "Schrodinger equation" 
idr^E,u — [tlE,u,^E.u] with tlic Hamiltoniau (TB] 

ilE,u - — { — +V2GfNA (1) 
Vu \2E ) 

+ —. — ^ / dE' / du' 

The matrix of neutrino mass-squares causes vacuum 
flavor oscillations and that of net charged lepton densi- 
ties — diag(rie— figj "■/^—'T'/i, ?^t~?^t) adds the Wolfen- 
stein matter effect. The third term provides neutrino- 
neutrino refraction and is analogous to matter except for 
Pantaleone's off-diagonal elements and except that in the 
SN context neutrinos are not isotropic. A neutrino ra- 
dial velocity at radius r is Vu = (1 — u E? /r'^ Y/'^ . The 
factor 1 — v^v^i arises from the current-current nature 
of the weak interaction and causes multi-angle effects. 
Moreover, Vu appears in the denominator because we fol- 
low the flavor evolution projected on the radial direction, 
causing the multi-angle matter effect |18j . 

Up to the MSW region, the matter effect is so large 
that ^E,u is very nearly diagonal in the weak-interaction 
basis, the usual approximation made in SN neutrino 
transport. Neutrinos remain stuck in flavor eigenstates 
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unless the off-diagonal ^e,u elements start growing by 
the self-induced instability. To find the latter we linearize 
the EoM in the small off-diagonal amplitudes. 

Stability condition. — We study the instability driven 
by the atmospheric Am^ and the mixing angle ^13 ^ 1, 
we work in the two-flavor limit, and switch to the lo — 
Am'^ /2E variable. We write the flux matrices in the form 



Tr*„ 



(2) 



where F^'f^ are the flavor fluxes at the inner boundary 
radius R. The flux summed over all flavors, Tr^^j.u, is 
conserved in our free-streaming limit. The ly^ survival 
probability is ^[1 + s^j^uif)] in terms of the "swap fac- 
tor" — 1 < Si^^uir) < 1. The off-diagonal element Suj,u is 
complex and „ + |5'^,„P = 1. 

The small-amplitude limit means liStj^ul 



linear order Su 



■C 1 and to 
1. Assuming in addition a large dis- 
tance from the source so that 1 — 1, the evolution 
equation linearized in Si^^u and in u is |23j 

-fi / du' duj' {u + u') g^'u' Sui',u' ■ (3) 



Weak-interaction effects are encoded in the r-dependent 
parameters with dimension energy 



A = x/2GF[ne(r)-ng(r) 



2r2 



V2Gf [F^^R) - F„JR)] R^ 

2r2 



(4) 



The factor i?^/2r^ signifles that only the multi-angle im- 
pact of the ly-v and matter effects are relevant for the sta- 
bility analysis, not the densities themselves. Both A and 
fi depend on R, but so does the occupied u-range: phys- 
ical results do not depend on the choice of R. We choose 
R — 44.7 km such that the occupied u-range is 0-1. We 
normalize the v-v interaction strength fi to the Ve-i^x flux 
difference at R, i.e. du dug^^u = —1, but the only 
physically relevant quantity is n{r)guj,u- Our SN model 
provides fJ-{R) = 1.73 x 10^ km^^ and an "asymmetry" 
£ — J duj du gi^ u = 0.35. 

Writing solutions of the linear differential equation, 
Eq. ([3]), in the form S^^u = e^'^'' with complex 

frequency ~ j + in and eigenvector Qui,u leads to the 
eigenvalue equation |23) . 



{lu+uX-Q) ( 



du duj' (u+u) gc^ 



(5) 



where A = A -I- e/i. The solution has to be of the form 
Qlu,u — {A + Bu)/{uj + uX~il). Solutions exist if — 
Ii ± \/Iol2, where /„ — ^ du) du g^ ^u'^ / {'^ + wA — fi). 
The system is stable if all Vl are purely real. A possible 
imaginary part, k, is the exponential growth rate. 
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FIG. 3: Growth rate k for our SN model as a function of fj, 
for various A values as indicated. 



Application to our SN model. — Ignoring the effect of 
matter (A = 0), we show k{^) for our 280 ms SN model 
in Fig. [3] The system is essentially stable above /i of a 
few 100 km~^. It is noteworthy that n is of the same 
order as a typical u of the g^^^u distribution, in our case a 
few km~^. We also show for A = 10^ and 10"^ km^^ 
and observe a shift to larger /i- values [23] . 

In Fig.|4]we show contours of k in the (/x. A) plane. For 
large /x and A values, the system is unstable for X ^ fi 
[23J. In other words, if the local neutrino number density 
is much smaller or much larger than the local electron 
density, the system is stable. 

We also show the locus of [/i(r), X{r)] along the radial 
direction. Since ^{r) cx r~^, the red solid line in Fig.|4]is 
essentially the SN density proflle. The step-like feature 
is the shock wave where the matter density drops by an 
order of magnitude. Without matter (A = 0), neutrinos 
would enter the instability strip at /i ~ 100, correspond- 
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FIG. 4: Contours for the growth rate a in km^^. Also shown 
is the profile for our SN model. The vertical axis essentially 
denotes the density, the horizontal axis the radius (fi oc r~*). 
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ing to r '-^ 150 km. We find similar results for other 
snapshots at times 150 and 400 ms postbounce, i.e., neu- 
trinos do not encounter a self-induced instability. 

Comparison with earlier work. — A similar accretion- 
phase model of the Basel group was used to study flavor 
stability by numerically solving the EoMs [in]- A mono- 
energetic neutrino distribution and isotropic emission at 
a neutrino sphere were assumed. For some snapshots, 
an instability occurred at a large radius. Applying our 
method to the same matter profile and schematic neu- 
trino distribution we find perfect agreement with Ref. |19j 
and even reproduce their onset radius for those cases 
where partial flavor conversion occurs |27j . It would be 
interesting to repeat our study with realistic Basel dis- 
tributions to see if partial flavor conversion is an artifact 
of their schematic energy and angle distributions. 

Conclusions. — We have performed a linearized flavor 
stability analysis of accretion-phase SN models and neu- 
trino fluxes with realistic energy and angle distributions. 
For these models, self-induced flavor conversions do not 
occur. One should apply this method to a broader class 
of models to see if this conclusion is generic. It also re- 
mains to extend a linearized analysis to cases without 
cylindrical symmetry of the angular distribution in view 
of Sawyer's concerns about a significant multi-angle in- 
stability [12j. In realistic 3D models, the neutrino dis- 
tribution is not cylindrically symmetric and even if this 
were the case, in principle even a small fluctuation could 
trigger a novel instability if it were to exist. 

Recent evidence suggests that the neutrino mixing an- 
gle 6*13 is not very small [2D1, implying that the MSW 
conversion in the SN is adiabatic. One can then dis- 
tinguish the neutrino mass hierarchy by Earth matter 
effects [3| or the early rise time pi , but only if collective 
oscillations do not swap flavors before the MSW region. 
(In the presence of the collective flavor swap, the hier- 
archy could be distinguished for an extremely small 9i3 
where the MSW conversion is not adiabatic [23 •) The 
suppression of self-induced flavor conversion during the 
accretion phase, if generic, is good news for the possibil- 
ity of determining the mass hierarchy with SN neutrinos 
if a "large" value for is experimentally conflrmed. 
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